Phase textures induced by dc current pairbreaking in multilayer structures and 

two-gap superconductors. 
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We predict current-induced formation of equilibrium phase textures for a multicomponent super- 
conducting order parameter. Using the two-component Ginzburg-Landau and Usadel equations, 
we show that, for weakly coupled co-moving superconducting condensates, dc current I first causes 
breakdown of the phase locked state at I > 7 c i followed by the formation of intrinsic phase textures 
well below the depairing current Id- These phase textures can manifest themselves in multilayer 
structures, atomic Bose condensate mixtures in optical lattices and two-gap superconductors, par- 
ticularly MgE>2, where they can result in oscillating and resistive switching effects. 
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PACS numbers: 74.20.De, 74. 20. Hi, 74.60.-w 

Interest in novel effects caused by multicomponent or- 
der parameters in heavy fermion and organic supercon- 
ductors has been recently amplified by the discovery of 
two-gap superconductivity in MgB2 , which exhibits an 
anomalous increase of the upper critical field by impuri- 
ties |3j and intrinsic Josephson effect between two weakly 
coupled order parameters = Aie 1 and if? 2 = Age a 
in a and ir bands. Excitations of the interband phase dif- 
ference 8 — 81 — 62 can be either the phonon-like Legget 
modes or phase textures formed by 2ir interband phase 
slips. Phase textures and peculiar vortex properties |f| 
in two-gap superconductors are generic manifestations 
of the multicomponent nature of the order parameter, 
which have analogs in superfluid 3 He 4:] and models of 
color superconductivity in the particle physics 6]. 

Spontaneous formation of phase textures breaking the 
time reversal symmetry is inhibited by interband cou- 
pling, which locks the phases of *f>i = Aie iSl and if? 2 = 
A2e 2 . However, for weak interband coupling character- 
istic of MgB2, the order parameters if?i and if? 2 can be 
unlocked by electric fields. The resulting non- equilibrium 
charge imbalance generates interband phase textures, 
which do not carry magnetic flux and thus do not interact 
with weak magnetic fields and supercurrents 0. In this 
Letter we show that, for sufficiently strong superconduct- 
ing currents, formation of phase textures does not nec- 
essarily require any nonequilibrium conditions and can 
result from dc current pairbreaking, which decouples if?± 
and ^2 well below the global depairing threshold. Thus, 
the phase textures are indeed a rather generic feature of 
equilibrium current-carrying states, which are not spe- 
cific to two-gap superconductors, but can be realized in 
any system with at least two different co-moving super- 
conducting or Bose condensates, in either the coordinate 
or the momentum space. The examples range from a 
weakly coupled thin film bilayer (Fig. 1) to a mixture of 
two weakly coupled atomic Bose condensates in optical 
lattices 8]. The bilayer in Fig. 1 can be mapped onto a 
two-gap superconductor in which two electron bands cor- 



FIG. 1: Formation of phase textures in current-carrying bi- 
layer. Thick black line shows the interlayer Josephson contact 
and gray rectangles show current leads. 



respond to the different films, and the interband coupling 
corresponds to the interlayer Josephson energy. 

The mechanism of texture formation is particularly 
transparent in a bilayer, which also reveals peculiar non- 
stationary effects such as switching, parametric amplifi- 
cation, and flux flow oscillations controlled by the applied 
current in the absence of a dc magnetic field. The sheet 
supercurrent I = I x + I 2 in the bilayer is sustained by 
the phase gradients 6[ and 8' in films 1 and 2, where 
h = -cdxfaOy&iPkl, I 2 = -cd 2 o 6>i/87r 2 A?, A x and 
A2 are the London penetration depths, 0o is the flux 
quantum, c is the speed of light, and the vector potential 
A is negligible for small film thicknesses d\ <C Ai and 
e?2 <C A2. For small current I, the minimum energy cor- 
responds to the phase-locked state with the global phase 
gradient 0[ = 8' 2 = Q = -87r 2 (A 2 /di + A 2 /d 2 )//c0 o for 
which the interlayer Josephson energy vanishes. 

The situation changes at higher Q ~ I/C2, if the co- 
herence length £2 in film 2 is greater than £1 in film 1. 
Then the Ginzburg-Landau (GL) current pairbreaking 
suppresses the gap A 2 (I) = A 2 [l — (^Q) 2 ] 1 ^ 2 more than 
the other gap Ai(7) = Ai[l - {£,iQ) 2 } 1/2 , so the max- 
imum phase-locked current is limited by the depairing 
threshold Q = 1/\/3& in film 2. For Q > 1/V3&, Aim 
2 goes normal, forcing current redistribution, which re- 
stores superconductivity in film 2, but causes a gradient 
of the phase difference 6'{x) and the oscillating perpen- 
dicular Josephson current JcSinf? along the bilayer. As a 
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result, current loops appear in which currents flow par- 
allel to I in the "stronger" film 1 and antiparallel to / in 
the "weaker" film 2, as shown in Fig. 1. Such decoupling 
transition occurs if the gain in the condensation energy 
due to the current redistribution exceeds the loss in the 
Josephson energy for a single vortex. 

A theory of this transition can be developed using the 
GL free energy = J FcPr, 



F =Y J d m [-a m A 2 m + >fAi + 7m (VA m ) 2 ] 
fS 2 /87r + e, 7 (l -cos9), 
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where m labels films 1 and 2, a m = a m — Q^jm, 
Qm = V0 m + 27rA/0o, and B = V x A. The last term 
defines the sheet coupling energy ej = 4>qJ c /2t:c where 
J c is the Josephson critical current density. For weak 
coupling, e(x) varies over the length Lg much greater 
than £ m = (7 m /a m ) 1/2 . Thus, A 2 , = a m /(3 m to the zero 
accuracy in ej, and the thermodynamic potential G for 
slowly varying Q m (x) and fixed /, takes the form: 
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Here Jj_ is the current density injected perpendicular to 
the layers, and the sheet current I(Q) along the bilayer 
is determined from OG/dA = 0. Next, we extract the 
0-dependent energy Gg of phase textures, expressing Q m 
in terms of e and / in Eq. |(2J). For slow variations of 
e(x), a quadratic gradient expansion of Gg yields 

r t 2 

Gg = e,, [^(V0) 2 + l-cos6-r]e-(QVe)h]d 2 r, (3) 



where r\ = J±/J c . The phase length Lg(Q), the coupling 
parameter h(Q) and I(Q) depend parametrically on the 
background gauge- invariant phase gradient Q: 
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/ = -Q[ £l (l - £ 2 Q 2 ) + e 2 (l - £Q 2 )]^c/^. (6) 

Here t\ = diQa\/2(3\ and e 2 = c^fa^/ 2/32 are charac- 
teristic condensation energies in films 1 and 2 at I = 0. 
For weak Josephson coupling, ej <C min(e m /£ 2 „), the 
phase length Lg is much greater than £i and £2, except 
for special cases discussed below. From Eq. © we obtain 
the following dynamic equation for 9(r,t): 

T 2 e + r r e = L 2 g V 2 e -sine- div(Q/i) + r), (7) 



where r 2 = C4>o/2ircJ c , r r 



)/2ttcRJ c , and C and R 



are the sheet capacitance and quasiparticle ohmic resis- 
tance of the Josephson contact, respectively. 



Eq. resembles the free energy of a long Josephson 
contact in a magnetic field, but in our case the driving 
term (QV#)/i results from the pairbreaking asymmetry 
of the layers. Here h{Q) cx I 2 at small I, but h(Q) di- 
verges at the global depairing current Id = I(Qd) where 
Qd = ( e i + £ 2)/3(ei£ 2 + £2^2)- Yet, Id cannot be reached 
because the phase-locked state becomes unstable above 
the depairing threshold, Q > l/\/3£2 < Qd, in film 2. 

At Q « the gradient term in Eq. changes 

sign so Gg should be expanded in higher order spa- 
tial derivatives of 6, which add the stabilizing term 
£ 2 L 2 {V 2 9) 2 /2 into Eq. @ where L 2 g {Q) « L 2 (1-3£ 2 Q 2 ), 
L = 2(e 2 /e,/) 1 / 2 , and £ ~ £ 2 - In the critical region, 
Q ~ 1/^3^2, a small perturbation 89 = 9q cos kx changes 
the energy by SGg cx 1 - k 2 L 2 (3&Q 2 - 1) + £ 2 L 2 k 4 . 
Hence, SGg is minimum at k 2 m = (3£ 2 Q 2 — l)/2£ 2 , and 
SGg{k m ) < if 3£fQ 2 > l + 2£/L . The bilayer becomes 
unstable with respect to phase perturbations with the 
wave vector k m = {£Lq)~ x I 2 at Q > Q C 2 where 
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At the spinodal point Q = Q c2 the phase-locked state 
is absolutely unstable, but stable large-amplitude tex- 
tures become energetically favorable at a lower Q w Q c \ 
due to proliferation of interlayer 27r phase slips 9 = 
4 tan -1 exv(x/Lg) similar to the Josephson vortices at 
H > H cl . At Q = Qd <C Qc2i the energy of a sin- 
gle phase slip 8ejLg equals the gain in the condensation 
energy 2nejQ c ih(Q c i) in Eq. 0. Hence, 



(ei + e 2 )ej 
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For I > J c i = 8ircQ c i(ei + £2)/4>o, the minimum energy 
corresponds to the chain of phase slips spaced by a(J) 
for which sin#/2 = cn(x/pLg\p 2 ) 0. Here cn(x|p 2 ) is an 
elliptic function, and 



2LgpK(p 2 ), I 3 =I 3 cl E(p 2 )/p, 



(10) 



where K(p 2 ) and E(p 2 ) are the complete elliptic inte- 
grals defined by the parameter < p < 1. The pe- 
riod a(I) diverges logarithmically at / — > I cl , but then 
a » 0.5n 2 Lg(I c i/I) 3 decreases rapidly as / further in- 
creases. As a result, the current density in layer 2 remains 
close to J2d, while J\ increases up to J\d, giving the max- 
imum current Id(T) = d\Jdi{T) + d2Jd2{T) for decoupled 
layers. Formation of the phase textures at Q = Q c i is 
therefore a first order phase transition with a spinodal 
decomposition of the phase-locked state at Q > Q C 2- 

The above GL theory of current-induced decoupling 
of spatially separated condensates holds if T c \ and T C 2 
are not too different, for example, in a bilayer made 
of the same superconductor with different concentra- 
tions of nonmagnetic impurities. It also implies that 
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T < T C 2 < T c i so superconducting states in both layers 
are weakly coupled, unlike the case T > T C 2 for which 
A 2 is induced by proximity effect, and the phase-locked 
state persists for all I < I d 0- Now we turn to the 
condensate decoupling in the momentum space, focusing 
on interband phase textures in two-gap superconductors. 
We use here the quasiclassic equations of two-gap super- 
conductivity in the dirty limit Q, which also enable us 
to calculate the phase textures in bilayers for all T: 
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(11) 
(12) 

(13) 



cj>0 m 



Here f m (r,uj) and g m (r,uj) are the Usadel functions in 
the m-th band, \f m \ 2 +g 2 n = 1, w = nT(2n + 1), D^J 3 are 
the intraband diffusivities, m = 2 if m = 1, and m = 1 
if in — 2, T m jn are the interband scattering rates, II = 
V + 2niA/(f> , N m is the partial density of states, X m m' 
are the BCS coupling constants, and iViAi 2 = N 2 \2i- 
The indices 1 and 2 correspond to the a and 7r bands of 
MgB 2 for which A„ « 3A W7r » 8A CT7r , AT W ps 1.3iV CT [3. 

Eqs. ifTTjl - lfT^l) can be obtained by varying the free 
energy 3 = J Fd 3 r where 



F =\J2 VANrnX-l, +F 1+ F 2 + F i + ^, (14) 



F m = 2ttT £ N m [u(l - g m ) - Re(/* * m ) 



O)>0 



+D^(n«/ m n^/* + V a g m Vpg m )/A\. (15) 

Here i 7 ^ is the intraband energy, and P; = 
2^TRe^ >0 (Ar 1 r 12 + Ar 2 r 21 )( 5l . 92 + /f/ 2 - 1) is due to 
interband scattering. The first term in Eq. (|14l) contains 
the Josephson-like interband coupling energy — qcos^ 
where e< = A x A 2 A 12 A r i/' u; , and to = A X1 A 22 - Ai 2 A 2i . 

We derive the equation for the slowly varying (Lg ^> £) 
interband phase difference 9(r), neglecting weak inter- 
band scattering 0] an d expanding Eqs. (fTT|> - <fT~Sf) in 
powers of Q 2 and Q 4 , where the term tx Q 4 accounts for 
current pairbreaking in the lowest order in J. This calcu- 
lation gives the energy of the phase texture in the form 
of Eq. (3J with t] = and 



2 _ 2(f 1 ip 2 
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h = 2Q^ 2iPl ~^ 2) 
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where Q = 87r 2 A 2 J/c0 o , A = o /[327r 3 (^i + <^ 2 )] 1/2 is 



the London penetration depth 
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" tm+ 32^[ tanh ^T + ^T SeCh ^ ' (18) 



Here u„ 



7rTA^E w >o^ 2 /(^ 
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27rTA 2 l ^ >0 (w 2 + A 2 l )- 3 / 2 . Eqs. (13-0 can also be 
applied to bilayers, by replacing J — > /, ej — > 4>oJ c /2irc, 
and iV m — * d m N m . Eq. reduces to Eqs. (QJ-® 
near T c if Q^ 2 <C 1. However, Eq. (|16fl also describes the 
case (particularly relevant to MgB 2 ) for which and 
^2 become weakly coupled only at low T as intraband 
pairing causes the Cooper instability in band 2. The sim- 
pler GL theory gives the full dependencies of Lg and h 
on current, while in the Usadel approach only the main 
quadratic term in h(Q) can be obtained analytically for 
all T. In particular, for T = 0, Eq. O yields 
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(19) 
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Here the band decoupling current density J c \ is de- 
fined as before by 8Lg = 2nQ c ih(Q c i) where h — 
(tt/4 + 4/37r)L 2 (£ 2 2 - £ 2 )Q 3 , Q cl = Sir 2 A 2 J cl /cfo, U = 
(Dm/Am) 1 / 2 is the intraband coherence length, and fc = 
(1/3 + 7r 2 /16) _1//3 « 1.017. If expressed in terms of the 
sheet condensation energy densities e TO = N m A 2 n ^ 2 n /2, 
Eqs. (fT^f) and (2Dj l reduce Eqs. © and © to the accu- 
racy of coefficients ~ 1 . 

The results presented above indicate that the previous 
calculations of Jd 0,0] for two-gap superconductors in 
which both bands were assumed to be phase-locked are 
only valid at higher temperatures T > T 2 ~ T c exp(l/A — 
I/A22) for which A 2 is induced by interband coupling. 
Here A = [A + + (A 2 _+4Ai 2 A 2 i) 1/2 ]/2u; and A± = An±A 2 2- 
To estimate T 2 for MgB 2 , we take T^=40K, An = 0.8, 
A 22 = 0.3, A12 = 0.12, A 2 i = 0.09 for which T 2 ~ 
0.12T C ~ 5K. This qualitative interpretation is consistent 
with the Usadel calculations and the strong coupling 
Eliashberg theory 0], which predict a two-hump J(Q) at 
low T, which turns into the conventional dome-like J(Q) 
at higher T. For T > T2, current pairbreaking thus occurs 
in the phase-locked state, for which J 4 = max(J(Q)). 

A two-hump J(Q) at low T means that four possi- 
ble phase gradients can provide the same current density 
J = J(Q). Since only states with d\ J\/dQ > are stable, 
the two-hump J(Q) would indicate formation of a strati- 
fied flow comprised of parallel channels with two different 
phase gradients, similar to the Gunn instability in semi- 
conductors. However, this stratification is preceded by 
interband decoupling, since for A12 <C A22, the band 2 
cannot sustain the same Q as the band 1. Thus, inter- 
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FIG. 2: Flux flow resistance R(I) for the moving phase texture 
as a function of the current I along the bilayer. 



band current redistribution provided by the phase tex- 
tures increases the GL depairing current density Jd(T) 
for which Qi = l/v3£i and Q2 = as compared 

to its phase-locked value Jd = maxJ((3): 



Jd 



Jd 



87tc(qi7i//3i + a 2 72//?2) 3/2 
3^ (7i//3i 



3\/30o 



0\ 



(21) 
(22) 



where a m , f3 m and 7 m are the 2-gap GL expansion coeffi- 
cients. The enhancement of Jd is most pronounced in the 
case of a clean a band (large 71 0) at low T < Th where 
Jd{T) could be measured by a pulse technique [LJ. 

Phase textures in bilayers manifest themselves in dc 
transport if both parallel and perpendicular currents / 
and J_l are applied, as shown in Fig. 1. For / < I c \, 
no voltage across the bilayer occurs, but for / > I cl , 
the oscillating voltage V(x,t) = — (cj) v/ircpLg)dn{(x — 
vt)/pLg\p 2 ] and the ohmic average voltage V = Rfl± 
appear due to the phase slip structure moving with the 
velocity v(I, I±). Here the resistance Rf is similar to the 
flux flow resistance of a long Josephson contact Q 



R f = Tr 2 R/4K(p 2 )E(p 2 ), 



(23) 



where R is the resistance of the interlayer contact, and 
p(I) is defined by Eqs. i|l(J|) . The dependence Rf on the 
longitudinal current 1 shown in Fig. 2 describes switch- 
ing between superconducting and resistive states across 
the bilayer. Other effects include a parametric resonance 
caused by superimposed ac currents I±(t) and I(t), since 
I(t) modulates the parameters in Eq. (0. Therefore, the 
geometry in Fig. 1 can provide switching and Josephson 
flux flow oscillator in current-operated devices. 

Phase textures in two-gap superconductors can move 
due to interband currents produced, for example, by 
nonequilibrium charge imbalance ■ One could also ex- 
pect a kink in the nonlinear low-frequency (w <C A 2 ) 
rf surface resistance as the field amplitude Hq exceeds 
the onset of the interband phase slip formation Hg = 



47rJ cl A/c. To estimate Hg for MgB2, we assume that 
the interband breakdown occurs as the screening cur- 
rent density J(O) ~ cHo/AttA on the surface produces 
the phase gradient Q — 8n 2 A 2 J(0)/c<^o exceeding the 
GL depairing limit Q C 2 = 1/\/3£t in n band. Hence, 
He = <Ao/2v3A^ 7r ~ U c ^ a j^, where H c is the ther- 
modynamic critical field, and the ratio £, a /£,-n can be 
strongly affected by impurities ||. Taking £ CT /£7r « 0.3 
for MgB 2 single crystals [3] and H c (0) ~ 0.3T, we find 
Hg(0) ~ 0.1T. Such fields cause breakdown of the lin- 
ear London electrodynamics, affect properties of vortex 
lattice, penetration vortices through surface barrier, etc. 

In conclusion, two weakly coupled co-moving supercon- 
ducting condensates can undergo a first order phase tran- 
sition into a phase textured state well below the global 
depairing current. Such textures controlled by current 
result in resistive switching and oscillating effects. 

This work was supported by US DOE Office of Science 
under contract No. W31-109-ENG-38. 

Note added. After this paper was submitted, a phase 
textured state with two different winding numbers in 
weakly coupled Al rings has been observed 
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